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We concern ourselves with the self-induced charge currents in electromagnetic materials and some 
related topics on its effects in the present paper. The contribution of self-induced charge currents of 
metamaterial media to photon effective rest mass is briefly discussed. We concentrate primarily on 
the torque, which is caused by the interaction of self-induced charge currents in dilute plasma with 
interstellar magnetic fields, acting on the torsion balance in the torsion balance experiment, ft is 
shown by our evaluation that the muons and alpha-particles in secondary cosmic rays will contribute 
an effective rest mass about 10 -54 Kg to the photon, which is compared to the newly obtained upper 
limit on photon rest mass in Luo's rotating torsion balance experiment. 

I. PHOTON EFFECTIVE REST MASS DUE TO SELF-INDUCED CHARGE CURRENT 

In is well known that in some simple electromagnetic media such as electron plasma, superconducting media and 
Lorentz dispersive materials (inside which the time-harmonic electromagnetic wave is propagating), the self- induced 
charge current density J [1] is proportional to the magnetic vector potential A. Thus the interaction term /j,qJ ■ A in 
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the Lagrangian density is therefore transformed into an effective mass term — \ m °^± A 2 of electromagnetic fields [2], 
which is analogous to the rest mass term in the Londons' electromagnetics for superconductivity, Ginzberg-Landau 
superconductivity theory and Higgs mechanism. Note that in these media, the magnetic permeability 1 /x = 1, so it 
is easy for us to obtain J (and hence the effective rest mass of photon) by solving the equation of motion of charged 
particles acted upon by electromagnetic waves. However, for some artificial composite electromagnetic metamaterials, 
which are described by, e.g., the two time derivative Lorentz material (2TDLM) model (in which both the electric 
permittivity and the magnetic permeability are of frequency dependence and therefore of complicated form) and some 
uniaxially (or biaxially) anisotropic media with permittivity and permeability being tensors, the above formulation 
may be not helpful in obtaining the effective rest mass of electromagnetic fields. In these cases, however, we have 
shown that the following formula [2] 

^ = (l-nV, (1-1) 

which arises from the Einstein-de Broglie relation, can be applicable to the above problem. In the expression (1.1), 
n, fi and c denote the optical refractive index, Planck's constant and speed of light in a vacuum, respectively. 

The two time derivative Lorentz material (2TDLM) model first suggested by Ziolkowski and Auzanneau [3,4] is 
a generalization of the standard Lorentz material model. Ziolkowski has shown that this type of 2TDLM medium 
can be designed so that it allows communication signals to propagate in the medium at speeds exceeding the speed 
of light in a vacuum without violating causality [5]. For these artificial metamaterials, the model encompasses the 
permittivity and permeability material responses experimentally obtained [6]. The x-directed polarization and y- 
directed magnetization fields in such materials would have the following forms [5] 

d 2 „ . „ d „ . a _ / 9 „ „ . „ d „ . „d 2 



^My + T^-My + ^My = ^Uy + U pX f ^Hy + ^ ^ Hy , (1.2) 



*This paper is concerned with the self-induced charge currents and some related physically interesting topics. It will be 
submitted nowhere else for the publication, just uploaded at the e-print archives. 
^E-mail address: jqshen@coer.zju.edu.cn 

1 Note that here for the magnetic properties in superconductors, we adopt the current viewpoint rather than the magnetic- 
charge viewpoint (where the permeability can be viewed as /i = 0). 
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where Xa' m j Xp™ an d Xj' m represent, respectively, the coupling of the electric (magnetic) field and its first and second 
time derivatives to the local electric (magnetic) dipole moments. ui p , T c , r m and uq can be viewed as the plasma 
frequency, damping frequency and resonance frequency of the electric (magnetic) dipole oscillators, respectively. Thus 
the frequency-domain electric and magnetic susceptibilities are given 

X — n i X — o i l -L -O 1 

-w 2 + iuT* + (w§) 2 -uj 2 + iuT™ + (wJJ 1 ) 2 

So, the refractive index squared in the 2TDLM model reads 

n 2 (w) = 1 + X° (w) + X m (w) + X° (w) X m (w) • (1.4) 
Thus it follows from (1.1) that the photon effective rest mass squared in 2TDLM model is of the form 

= ^ [x° M + x m M + x c M x m M] • (1.5) 

As an illustrative example, we apply Eq.(1.5) to the case of left-handed media, which is a kind of artificial compos- 
ite metamaterials with negative refractive index in the microwave frequency region (10 9 Hz). More recently, this 
medium captured considerable attention in various fields such as condensed matter physics, materials science, applied 
electromagnetics and optics [6-11]. The most characteristic features of left-handed media are: (i) both the electric 
permittivity e and the magnetic permeability n are negative; (ii) the Poynting vector and wave vector of electromag- 
netic wave propagating inside it would be antiparallcl, i.e., the wave vector k, the electric field E and the magnetic 
field H form a left-handed system; (iii) a number of peculiar electromagnetic and optical properties such as the reversal 
of both the Doppler shift and the Cherenkov radiation, anomalous refraction, modified spontaneous emission rates 
and even reversals of radiation pressure to radiation tension [7,8] arise. All these dramatically different propagation 
characteristics stem from the sign change of the group velocity. In what follows we calculate the photon effective rest 
mass squared by making use of Eq.(1.5) and the expressions 

uj 2 Flu 2 
e(u) = 1- . p . , , mH = 1--2 g -f (1.6) 

w {uj + ii) v ' to 2 - col + icoT 

for the effective dielectric parameter and magnetic permeability, where the plasma frequency uj p and the magnetic 
resonance frequency u>o are in the GHz range {e.g., lu p = 10.0 GHz, loq = 4.0 GHz [12]), the damping parameters 
7 = 0.03w p , T = 0.03^0- The parameter F may often be chosen < F < 1, for example, F — 0.56 [12]. By using 
n 2 (uj) — e (lu) fj, (u), we obtain [13] 



It follows that the mass squared m 2 ff is obtained 
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off 



(1.8) 



Note that for the electron-plasma medium, where the electromagnetic parameters 7 = T = 0, F = 0, one can arrive 
at the following familiar formula 

m 2 cS = (1.9) 

with the plasma frequency squared w 2 = (here e, m c and N respectively denote the electron charge, mass and 
electron number density in this electron plasma). 

It should be emphasized that there is a difference in the definition of photon effective rest mass between two types 
of experimental schemes of testing photon mass 2 . In the discussion of wave dispersion and potential variations of the 



2 These two types of experimental schemes are as follows: one is based on the Hamiltonian density of electrodynamics and the 
other based on Maxwellian equations. 
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speed of light with frequency, which are based on the Maxwellian equations (rather than the Hamiltonian or Lagrangian 
density of interacting electromagnetic system), the effective rest mass squared of photons is expressed by Eq.(1.9). 
Historically, this problem has been discussed by many investigators, e.g., Feinberg who considered the possibility of 
the variation of the speed of light with frequency in the sharply defined optical and radio pulses from pulsars and took 
into account in more detail the observed variation of arrival time with frequency for the radio waves attributed to 
the interaction with interstellar electrons [14]. However, for those experimental schemes of testing photon rest mass, 
which are based on the Hamiltonian density (rather than the Maxwellian equations), the effective rest mass of photons 
due to media dispersion (i.e., the interaction of wave with charged particles in environmental dilute plasma such as 
the secondary cosmic rays) is two times that expressed in Eq.(1.9). The torsion balance experiment is just this type 
of schemes, which will be considered in what follows. 



II. EFFECTS OF SELF-INDUCED CHARGE CURRENTS IN THE TORSION BALANCE EXPERIMENT 



In this section, we consider the potential effects of self-induced charge currents arising in the torsion balance 
experiments [15,16]. The Lagrangian density of electrodynamics reads 



l -F F 1 

"4 M M!y 2 



F^F^ - -n^A^ + mo (2.1) 



where p 2 = (-^-) 2 with m 7 , % and c being the photon rest mass, Planck's constant and speed of light in a free 
vacuum, respectively, and Mo denotes the magnetic permeability in a vacuum. It follows from (2.1) that the canonical 
momentum density is written 

dC 

^ = ^T' ^) = -E(z). (2.2) 

In the electron plasma, due to the conservation law of the canonical momentum density, i.e., ^(to c v + eA) = 0, we 
have v = — ^-A + C with C being a constant velocity. By using the formula H = — E • A — C and electric current 
density J = New = — ^^A + NeC, one can arrive at (in SI) 



«-r- 

2/x 



^ + B^A 2 + 2^A2 + ^(v.E-^ 

c 1 m c p z c l \ eo 



JoAo-iVeC-A, (2.3) 



where Jo = p and A = </> are respectively the electric charge density and electric scalar potential, and e represents 

the electric permittivity in a vacuum and c = -y==- Note that here -^s • E — results from both the mass 
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term ^Aq and the Gauss's law V • E = -^A + 

It follows from (2.3) that the total effective rest mass squared of electromagnetic fields in electron plasma is given 
as follows 

= + (2-4) 
iti e 

It is worthwhile to point out that according to the Ampere-Maxwell-Proca equation V x B = Mo^oJjE + p J — M7A 
and the consequent VxB = Moeo J^E- [fi 2 , + ^f ]A + p Q NeC, the p, 2 ot should be p 2 ot = p 2 + rather than that 

in (2.4). This minor difference between these two p 2 ot results from the derivative procedure applied to the Lagrangian 
density by using the Euler-Lagrange equation. Although it seems from the Ampere-Maxwell-Proca equation that 
p 2 ot = p 2 + t *"J^ e , the factor-2 in (2.4) cannot be ignored in calculating the torque acting on the torsion balance 
due to the photon effective rest mass, which will be confirmed in what follows. In the classical electromagnetics, the 
following familiar formulae are given 

V x B = /i J, m d = nr 2 I, I = J • S, f=m d xB (2.5) 

and 

VxA = B, a d = nr 2 — , $ = B ■ S, (2.6) 
Mo 
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where V x B = /in. J holds when the electric field strength E, electric current density J are vanishing in the Ampere- 
Maxwell equation. In Eq.(2.5), m<j = irr 2 I means that a current loop of radius r carrying current I gives rise to a 
magnetic dipole moment m^; S in I = J • S denotes the area vector through which the electric current density vector 
J penetrates; r = x B means that an electric-current loop with a magnetic dipole moment m^, immersed in a 
magnetic field B, experiences a torque f = x B. Eqs.(2.6) means in the torsion balance experiment, a toroid coil 
contains a loop of magnetic flux $ which acts as a dipole source ad of magnetic vector potential via Vx A = B with 
the magnetic field B within the toroid as the source term. 

It follows from (2.3) that in the Hamiltonian density Tt of electromagnetic system, there are ^B 2 and 

\ (ji 2 + 2 Mn T ^ re ^ A 2 (i.e., | (/^totA) 2 ), which provide us with some useful insights into calculating the torque produced 

by the interaction of magnetic dipole vector potential moment with the ambient cosmic magnetic vector potential 
acting upon the toroid in the torsion balance experiment. For convenience, the Equations (2.6) may be rewritten as 
follows 

V x (MtotA) = no , (MotCLd = 7rr , = • S. (2.7) 

V Mo / V Mo / Mo V Mo / 

In order to calculate the torque acting on the torsion balance, we compare the equations and expressions in (2.7) 
with those in (2.5) as follows: 





the familiar case 


the case due to photon mass 


inH 


IB 2 


\ (/i to tA) 2 


Equations 


V x B = fi J 




dipole moment 


rrid = nr 2 I 


» tot a d = nr 2 (^) 


flux 


I = J • S 




torque 


r = nid x B 


f =? 



*Table 1 



Thus it follows from the analogies in the above table that, ? = (/itot a d) x (MtotA), i.e., the torque f = /i 2 ot &d x A. 
Here the torque f arising from the interaction between the dipole field of magnetic potentials (resulting from the 
electric current I in the toroid windings) with the ambient vector potential A acts on the torsion balance. It should 
be noted again that here /i 2 ot = /i 2 + rather than /i 2 ot = [i 2 + 

Note that there exists a term fi NeC in the Ampere-Maxwell-Proca equation V x B = /i eo J^E — [/i 2 + A "^ e ]A + 
HoNeC. Does this term has influences on the torsion balance? In what follows we will discuss this problem. 

In the similar fashion, we consider the familiar Hamiltonian density M • B, which describes the interaction of 
the magnetic moments with the magnetic field B. According to the Ampere-Maxwell-Proca equation with ^ £o JjE 

and tlaNe A being ignored, we have V x B = u NeC, which can be rewritten as V x B = (uoMtot) ). In the 

meanwhile, H — —NeC ■ A can also be rewritten as H = — (^^f^j ' (MtotA). 

It is known that in the sufficiently large space filled with homogeneous magnetic moments with volume density 
(magnetization) being M, the interior magnetic field strength B produced by homogeneously distributed magnetic 
moments is just /ioM 3 , i.e., V x A = ^i M or V x (yU to tA) = (/ioMtot) M. 



3 This may be considered in two ways: (i) if the direction of M is assumed to be parallel to the third component of Cartesian 
coordinate system, then the potentially non-vanishing components of B and H are also the third ones. It follows from the 
Ampere's law V x H = Jf ree with Jf ree = that H3 is constant, and Hz, can be viewed as zero. So, according to the expression 
B = /ito(H + M), we obtain B = /UoM; (ii) In accordance with Ampere's circuital law, the magnetic induction B inside a long 
solenoid with JV circular loops carrying a current / in the region remote from the ends is axial, uniform and equal to /io times 
the number of Ampere-turns per meter nl, i.e., B — /ion/, where n = y with I being the solenoid length. Let S denote the 
cross-section area of this long solenoid. The volume density of magnetic moments IS resulting from the carried current / is 
M = UJilllL with the solenoid volume V = SI. This, therefore, means that the volume density of magnetic moments just equals 
the number of Ampere-turns per meter nl, i.e., M = nl. So, we have B = /iqM. 
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Thus the comparisons between "the familiar case" and the case due to the constant C are illustrated in the following 
table, 





the familiar case 


the case due to constant C 


in H 


±B 2 


i (fl tot A) 2 


interaction Hamiltonian density 


h = -m-b 




Equations 


V x (MtotA) = Oo^tot) M 


VxB = ( M0 Mtot)(^) 


torque 


f = M x B 





*Table 2 



It follows that ? = ^-^j^J x (/i to tA), i.e., r = iVeC x A, which dose not act on the toroid in the torsion balance 
experiment. 

Since we have discussed the related preliminary preparation for the effects of self- induced charge currents, we think 
the order of magnitude of its effects in experiments {e.g., the extra torque due to the photon effective mass acting on 
the toroid in the torsion balance experiment) deserves evaluation. 



III. DISCUSSION AND REMARKS: SELF-INDUCED CHARGE CURRENTS IN DILUTE PLASMA 

Here we discuss the effective rest mass of photon resulting from the self-induced charge currents in the torsion 
balance experiments [15,16]. The self-induced charge current arises mainly from the two sources around us: (i) 
muon (/x) component and alpha-particles in cosmic rays; (ii) decay daughter (alpha-particles) of radioactive radon gas 
(|g 2 Rn, ||°Rn) in room environment. 

It is known that at the sea level, the current density of muon component in secondary cosmic rays is about 
1 x 10~ 2 cm~ 2 • s _1 . Assuming that the muon velocity approaches speed of light, the volume density of muon can be 

derived and the result is = 0.3 x 1CP 6 m~ 3 . So, according to Eq.(2.4) {i.e., m c g = f \J ia 2 f e ^ ) , electromagnetic 

wave with wavelength 4 A 3> 100 m at the sea level acquires an effective rest mass about 0.6 x 10 -53 Kg. Note that the 
above-mentioned muon volume density (0.3 x 10~ 6 m~ 3 ) is the datum at the sea level. In the Luo's rotating torsion 
balance experiment whose total apparatus is located in their cave laboratory, on which the least thickness of the cover 
is more than 40 m [15], it is reasonably believed that the muon density of cosmic rays in the vacuum chamber of Luo's 
experiment may be one or two orders of magnitude 5 less than that at the sea level, and therefore the effective rest 
mass acquired by photons in the vacuum chamber of Luo's torsion balance experiment may be about (0.6 ~ 2) x 10~ 54 
Kg. 

It is believed that the air molecules in the vacuum chamber of torsion balance experiment cannot be ionized by the 
alpha-particles of cosmic rays, because the mean free path of alpha-particle moving at about 10 7 m • s _1 in the dilute 
air with the pressure being only 10~ 2 Pa [15] is too large (more than 10 6 m) 6 . So, the air medium has nothing to 
contribute to the effective rest mass of photons. As far as the alpha-particle itself is concerned, it follows, from the 
fact that at the sea level the alpha-particle number density in secondary cosmic rays is about one order of magnitude 
less than that of muon, that the alpha-particle in cosmic rays itself contributes about 10~ 54 Kg to photon rest mass. 



4 Since the mean distance of two muons in secondary cosmic rays at sea level is about 100 m, Lorentz's mean-field formulation 
(Lorentz's dispersion theory, Lorentz's electron theory, 1895) is not applicable to the electromagnetic wave with wavelength 
A <C 100 m. So, the above effective mass formula is no longer valid for this case {i.e., A <C 100 m). In fact, the electromagnetic 
wave with wavelength A <C 100 m dose not acquire this effective rest mass. But the ambient cosmic magnetic vector potentials 
(interstellar magnetic fields) with low (or zero) frequencies will truly acquire this effective rest mass. 

5 Readers may be referred to the following data of muon current density in underground cosmic rays: the muon current densities 
are 10 _4 cm~ 2 • s _1 and 10 _6 cm -2 • s _1 at the equivalent water depth of 100 m and 1000 m under the ground, respectively. 

6 According to some handbooks, the free paths of alpha-particles in air (with pressure 10 5 Pa) are 2.5, 3.5, 4.6, 5.9, 7.4, 8.9, 
10.6 cm, corresponding respectively to the energy 4.0, 5.0, 6.0, 7.0, 8.0, 9.0, 10.0 Mev. So, in the vacuum chamber with pressure 
10~ 2 Pa in Luo's experiment, the free path of alpha-particles (with energy much more than 10 Mev, even perhaps more than 
100 Mev) in secondary cosmic rays is more than 10 7 times several cm, which means that the air molecules in the low-pressure 
vacuum chamber cannot be easily ionized by the alpha-particles of cosmic rays. 
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Radon (sgRn) possesses the a-radioactivity. Generally speaking, the Rn density in room environment is 10 -8 m~ 3 or 
so, which may be the same order of magnitude as particle density in secondary cosmic rays , and the decay daughter 
(alpha-particles) will therefore give rise to a photon effective mass of about 10~ 54 Kg. 

Thus the effective rest mass induced by muons and alpha-particles in secondary cosmic rays and radon gas can be 
compared to the new upper limit (1.2 x 10~ 54 Kg) obtained by Luo et al. in their recent rotating torsion balance 
experiment [15]. So, as far as the contribution of ions of cosmic rays to the photon rest mass is concerned, this upper 
limit (1.2 x 1CP 54 Kg) becomes just a critical value. We think Luo's experimentally obtained upper limit is of much 
interest and therefore deserves further experimental investigation to get a more precise upper limit on photon mass. 

IV. CONCLUSION 

In this paper, 

2 4 

i) it is verified that the expression (1.1), i.e., m ° f ^ c = (1 — n 2 )uj 2 can be applied to calculating the effective rest mass 
of photons in the artificial composite electromagnetic metamaterials with complicated permittivity and permeability 

ii) it is shown that the extra torque 2 Mo ^ e x A by the self-induced charge currents will also acts upon the toroid 
in the torsion balance experiment whose apparatus is immersed in the dilute plasma such as the cosmic rays and 
radon gas. Here the photon effective rest mass squared is fi 2 s = 2 t±0 ^ e rather than ^ ff = ^"^ e . 

The following difference between these two [i 2 s should be pointed out: the theoretical mechanism of torsion balance 
experiment lies in the Lagrangian or Hamiltonian (2.3) of interacting electromagnetic system, so, the effective rest 
mass squared of electromagnetic fields is /i 2 g = 2 Mo t ^ fe . However, for those experimental schemes [14,17,18] of photon 
rest mass stemming from the Ampere-Maxwell-Proca equation V x B = (U e o^E— [/j, 2 + Mo ^ e ]A + noNeC, the effec- 
tive rest mass squared of electromagnetic fields is n 2 s = rather than /j, 2 s = 2 ^ oJVe . The so-called experimental 
realizations stemming from the Ampere-Maxwell-Proca equation are as follows: pulsar test of a variation of the speed 
of light with frequency [14]; geomagnetic limit on photon mass based on the analysis of satellite measurements of the 
Earth's field [17]; experimental test of Ampere's law at low temperature [18], etc.. In all these experimental schemes, 
the effective mass squared is fj, 2 s — A "^ e . 

iii) the self- induced charge currents in dilute plasma presented above (i.e., muons and alpha-particles in secondary 
cosmic rays and radon gas) will give arise to an effective rest mass of about 10 -54 Kg, which can be compared to the 
new upper limit on photon rest mass obtained in Luo's rotating torsion balance experiment. 
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